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Progress Report 


The utility of the Mars Planetary Boundary Layer Model (MPBL) for 
calculations in support of the Mars 94 balloon mission has been 
substantially enhanced by the introduction of a balloon equation of motion 
into the model. Both vertical and horizontal excursions of the balloon are 
calculated along with its volume, temperature, and pressure. The 
simulations reproduce the expected 5-minute vertical oscillations of a 
constant density balloon at altitude on Mars. Figures 1-3 show the results 
of these calculations for the nominal target location of the balloon 
(latitude 50 degrees; solar longitude 170 degrees; thermal inertia 360 SI; 
visible optical depth 0.45; infrared optical depth 0.30; surface albedo 
0.20; surface pressure GOOmb) referred to as Moyen 50m by the French and 
Russians. 

Papers giving full details of the MPBL and its coupling to the Martian 
regolith for water vapor transport calculations are now In preparation. A 
paper entitled "A coupled subsurface atmosphere boundary layer model of H 2 0 
on Mars" by A. P. Zent, R. M. Haberle, H. Houben, and B. M. Jakosky was 
presented at the 22nd Lunar and Planetary Science Conference in March. The 
abstract is attached as Appendix A. 

A nonlinear balanced model (renamed the Balance System Model to emphasize 
the fact that a full dynamically-consistent set of equations is solved, 
rather than just the balance equation) has been developed for the Martian 
atmosphere. This model takes temperature data (such as are to be expected 
from satellite missions) and computes dynamically consistent rotational 
wind fields, dlabatic heating rates, and diabatic circulations. Results of 
this model have been compared with those from the full primitive equation 
model calculations of the Mars Climate Model (MCM) under both clear and 
dusty conditions. The accompanying Figures 4-7 illustrate the Balance 
System Model's ability to diagnose winds using only information about the 
mass field. 

The Balance System Model has been used to initialize a primitive equation 
model for simulations of the earth's stratosphere at the time of the El 
Chichon eruption in 1982. (The available satellite observations include 
only the temperatures and thus closely resemble the expected Martian 
observations.) As is the case with many such model simulations, a cold 
pole develops over tlme--probably due to the inadequate representation of 
small-scale gravity waves. Therefore, the Balance System Model is used as 
an assimilation model to update the temperature and wind fields at frequent 
intervals. The results for May and June (the third and fourth months of 
the simulation) are illustrated in Figures 8-15. 

Future Work 

Along with the completion of the papers now in preparation, work on the 
MPBL will emphasize those elements needed to further delineate the slope 
wind model. In particular, a self consistent pressure calculation, while 
computationally unwieldy, should allow a better forecast of slope winds 
(possibly even involving complex terrains in a 2-dlmensional model). This 
improved modeling, along with further study of the newly available 
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topography and figure data for Mars, should result in better simulations of 
the Viking Lander meteorological data. Of particular concern is the very 
strong vertical wind shear inferred from the Viking entry science. 


It is expected that further extensions of the current suite of Mars models 
to the other terrestrial planets will be undertaken in the near future. Of 
particular interest is the planet Venus, where recent spacecraft missions 
(Galileo and Magellan) have generated much excitement about that planet's 
atmosphere and surface-atmosphere interactions. 
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Figure Captions 


Figure 1. MPBL results for the nominal target point of the Mars 94 balloon 
mission. (See text for particulars.) Plotted as functions of time of day 
are balloon altitude, balloon and ambient temperatures, and balloon and 
ambient pressures. 

Figure 2. Further MPBL results using the balloon equation of motion. The 
vertical velocity and balloon overpressure plots show clear indications of 
the 5 minute vertical oscillation which is characteristic of a constant 
volume balloon under mars conditions. 

Figure 3. The horizontal displacement of the balloon as a function of time 
of day (symbols are plotted every hour) based on the above MPBL 
calculations. A significant fraction of the balloon’s travel takes place 
at night while it is "parked" at low altitude in spite of the frictional 
drag which must be considered under those circumstances. 

Figure 4. A comparison of winds generated by a 3-D primitive equation 
model of the Martian atmosphere (MCM) in a dust free simulation with the 
winds diagnosed by the Balance System Model using only the temperature 
fields of the MCM. a) The zonal mean temperature field of the MCM; b) the 
diagnosed diabatic heating rates (degrees per day) in the atmosphere; c) 
the MCM mean zonal wind; and d) the mean zonal wind diagnosed by the 
Balance System Model based only on the temperatures in (a). 

Figure 5. A continuation of the comparison between the MCM and the Balance 
System Model in a dust-free atmosphere, a) The mean meridional circulation 
of the MCM (the density-weighted streamlines shown are in units of 0.1 
megatons/second of material transport across a latitude circle); b) the 
mean meridional circulation diagnosed by the Balance System Model; c) 
streamlines of the horizontal flow for the MCM (units are m/s) at the 3mb 
pressure level; and d) the corresponding streamlines of the Balance System 
model . 

Figure 6. Another comparison between the MCM and Balance System Model as 
in Fig. 4, but for a dusty atmosphere with optical depth 1. 

Figure 7. As in Fig. 5, but for dust optical depth 1. 

Figure 8. A comparison between the observed average stratospheric 
temperatures in May 1982 (from the National Meteorological Center) and 3-D 
model calculations. Since the mean temperature field in the model Is 
updated using the Balance System Model as an assimilation routine, this 
figure illustrates the drift between model and observed temperatures over a 
period of 24 hours. 

Figure 9. a) The computed mean zonal wind in the stratosphere in May 1982. 
This quantity is not directly observed. The developing south polar night 
jet and the summer hemisphere easterly mesospheric jet are quite evident, 
b) The computed mean meridional wind in the stratosphere in May 1982. All 
winds are quite light except in the highly viscous sponge layer near the 
top of the model. Units are m/s. 
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Figure 10. A comparison between modeled and observed geopotential height 
fields in the stratosphere in May 1982 for zonal wavenumber one. The model 
generally shows larger amplitudes than the observations, but agreement 
below lOmb (the region of greatest interest for considerations of transport 
of the El Chichon cloud) is good. 

Figure 11. A comparison between modeled and observed geopotential height 
fields in the stratosphere in May 1982 for zonal wavenumber 2. There is 
some indication in the observations of instability of the polar night jet 
at high altitudes. Strong friction in the model at these heights damps any 
such instabilities. 

Figure 12. As in Fig. 8, for June 1982. 

Figure 13. As in Fig. 9, for June 1982. The polar night jet is continuing 

to accelerate, consistent with climatological data. 

Figure 14. As in Fig. 10, for June 1982. 

Figure 15. As in Fig. 11. for June 1982. In this plot there is the hint 

that the instability of the polar night jet is overcoming even the strong 
damping of the model sponge layer. 
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Appendix A 


A COUPLED SUBSURFACE ATMOSPHERE BOUNDARY LAYER MODEL 
OF H 2 0 ON MARS, A P. Zent, SETI Institute and NASA Ames Research Center, R. H. 
Haberle, NASA Ames Research Center, Moffett Field, Ca 94035, H. Houben. Space Physics Research 
Institute, Sunnyvale, Ca. 94087, B. M. Jakosky, LASP, University of Colorado, Boulder, Co. 80309. 


We have developed a coupled subsurface* atmosphere boundary layer model of H 2 0 exchange 
on Mars. Our objectives in constructing this model are: a) to identify the physical processes that 
control exchange of water between the atmosphere and regolith over a diurnal timescale on Mars, 
and to define physical parameters that can be used predictively; b) to understand the intensity as 
well as temporal and spatial variability of the H 2 0 flux through the martian surface; c) to develop 
a description of how exchange patterns vary throughout reasonable parametric space; and d) to 
advance our understanding of the martian annual H 2 0 cycle. 

Model: The atmospheric boundary layer model is taken from Haberle tt ai 51990, and is 
described in detail by those authors. It is a radiative-convective model which allows for the radiative 
effects of dust. The model is similar to that by Flasar and Goody (1976), however the model has 
an improved treatment of turbulence, and provides a self-consistent calculation of surface stresses, 
which help regulate the exchange of water between the surface and atmosphere. 

The subsurface model consists of a thermal component and a H 2 0 transport component. 
The heat conduction equation was solved in order to find subsurface temperatures. Viking data 
was used to constrain the thermal constants. Boundary conditions are found from the atmosphere 
model, which calculates the radiative and sensible heat fluxes. The surface temperature is then 
calculated assuming that the radiative, sensible, conductive, and latent heat fluxes sum to zero. The 
lower boundary condition is found by extending the solution region obtained to depths greater than 
the penetration depth of the annual thermal wave. The thermal gradient at the lower boundary is 
then assumed to be the quotient of the heat flux, calculated assuming a chondritic Mars, and the 
conductivity. 

There are four assumptions that underlie the model of H 2 0 transport through the regolith. 
The first is that vertical transport of H 2 0 through the regolith is according to Fick’s law. We 
further assume that water in the regolith must exist in one of three discrete states, either as vapor, 
adsorbate, or ice. 

<r = h + 0(7 >T) + * (1) 


Where a is the total density of water per cubic meter of regolith, / is the porosity of the 
regolith, 7 is the vapor density per cubic meter of gas , a is the adsorbed phase, and i is the ice 


phase. 

Some assumption must be made regarding the form of the expression 0(7, T) 
assumption is that 

«(7 ,T) = 
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( 2 ) 


Where 0 = 2.043x 10“*, and 6 is -2679.8 (Fanale and Cannon, 1971). P is the partial pressure of 
H 2 0 , leading to our final assumption, that the ideal gas law is obeyed by H 2 0 . 

The evolution of the population in a finite volume of regolith is controlled by 


da _ d<f> _ d dy 
c tt dz dz dz 


( 3 ) 


Where a is the total concentration of H 2 0 in the regolith, 7 is the vapor phase concetration, 
and <t> is the flux. The upper boundary condition is a flux, which is calculated by assuming that the 
vapor pressure of H 2 0 at the mathematical surface is controlled by adsorption. That permits the 
atmosphere layer model to calculate a gradient and resulting flux through the surface via stability 
functions from Monin-Obukhov similarity theory. The lower boundary condition for H 2 0 is zero 
flux. 

The model is initialized with no H 2 0 in the atmosphere and 2.1 kg m” 3 throughout the 
regolith. We carry out the calculation until the day to day variation of H 2 0 in the atmosphere 
column is less than 0.01%. Accordingly, the atmospheric H 2 0 column is stable over a timescale 
longer than the diurnal radiation pattern remains constant. We note that equilibration typically 
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requires 12 to 18 sols, providing the first estimate of the timescale for the atmosphere and regolith 
to equilibrate. 

Results: As our baseline case, we take the latitude and season that the Viking One lander 
touched down; albedos and thermal inertias are consistent with the Viking data. 

Beginning at midnight, the atmosphere is quite stable, the only turbulence is due to wind 
shear (Fig. 1). The regolith cools progressively throughout the night, increasing its adsorptive 
capacity (Eq. 2 and Fig. 2). The water vapor mixing ratio in the lower part of the atmosphere 
decreases as the night progresses because the cooling regolith is scavenging H 2 O ; due to the stability 
of the atmosphere, supply of water to the regolith is effectively diffusion limited. 

Sunrise occurs just before 0600 hours, and the soil begins to warm. Although the atmosphere 
is not yet convective, there is a thin layer of atmospheric mixing that brings additional water 
down to the surface, where the warming regolith continues to adsorb more water. Although physical 
adsorption is inhibited at higher temperatures, we predict continued adsorption after sunrise because 
the regolith is not in equilibrium with the atmosphere as a whole, but only with the strongly depleted 
lowest tens of meters. 

Even after the atmosphere becomes fully convective, about 0800 hours, there is still a short 
period where the atmosphere is supplying water to the regolith. The flux into the soil increases 
dramatically when the atmosphere begins to convect freely (Fig. 3). It isn't until about 1000 hours 
that the increasing temperature and adsorbed water content raise the pore pressure of water to 
atmospheric levels and the flux reverses. 

During the day the atmosphere is fully convective, and there is only a very shallow vertical 
gradient in atmospheric water. The regolith communicates during the day with the entire boundary 
layer, and approximate equilibrium exists. Even before the surface temperature peaks, the flux to 
the atmosphere begins to drop as the amount of near-surface HjO is depleted. Again in the late 
afternoon, there is a strong pulse of H 2 0 into the regolith, as the atmosphere is still convecting 
vigorously, but the regolith is now cooling rapidly. Finally, once the surface cools below the tem- 
perature of the atmosphere, about 1650 hours, convection stops, the atmosphere stabilizes, and the 
slow scavenging of HjO from the lowest tens of meters of the atmosphere begins again. 

We have examined the effects of albedo, thermal inertia, latitude, atmospheric optical depth, 
and pore size on our baseline model, and will report those results as time permits. 
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